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I present the preliminary results of a program to measure the star formation history in the halo 
of the Andromeda galaxy. Using the Advanced Camera for Surveys (ACS) on the Hubble Space 
Telescope, we obtained the deepest optical images of the sky to date, in a field on the southeast 
minor axis of Andromeda, 51' (11 kpc) from the nucleus. The resulting color-magnitude dia- 
gram (CMD) contains approximately 300,000 stars and extends more than 1.5 mag below the 
main sequence turnofi', with 50% completeness at = 30.7 mag. We interpret this CMD using 
comparisons to ACS observations of five Galactic globular clusters through the same filters, 
and through x'^-fitting to a finely-spaced grid of calibrated stellar population models. We find 
evidence for a major (~ 30%) intermediate-age (6-8 Gyr) metal-rich ([Fe/H]> —0.5) population 
in the Andromeda halo, along with a significant old metal-poor population akin to that in the 
Milky Way halo. The large spread in ages suggests that the Andromeda halo formed as a result 
of a more violent merging history than that in our own Milky Way. 



1. Introduction 

One of the primary quests of observational astronomy is understanding the formation 
history of galaxies. An impediment to this research is the relative paucity of galax- 
ies in the Local Group, which contains no giant ellipticals, and only two giant spi- 
rals - our own Milky Way and Andromeda. Fortunately, Andromeda (M31, NGC224) 
is well-situated for studying the formation of giant spiral halos, due to its proximity 
(770 kpc: 'Freedman fc Madore 199fl|l . small foreground reddening [Eb-v = 0.08 mag; 
ISchlegel' Finkbciner, fc Davis 1990| |, and low inclination [i « 12.5°; dc Vaucouleurs 195§^. 
Andromeda is similar to the Milky Way in many respects (Hubble type, absolute mag- 
nitude, mass, and size; |van den Bergh 1992| [Klypin, Zhao, k, Somerville 2002| ), but we 
have long known that its halo is more metal-rich than that of the Milky Way; the metallic- 
ity distribution in the Milky Way halo peaks near [Fe/H]Ri -T 8 dRyan fc Norris 1991]), 
while that in the Andromeda halo peaks near [Fe/H]« -0.5 rMould Kristia,n 19861 
[Holland, Fahlman, fc Richer 1996 Durrell, Harris, fc Pritchet 2001 ). Ahhough the Milky 
Way halo is dominated by old stars (e.g., f^anSenSe^^OOOt , the formation history of 
the Andromeda halo has been unknown until now. 

Physical processes possibly at work in forming spiral halos include rapid dissipative 
collapse in the early u niverse (|Eggen, Lynden-Bell, fc Sandage 1962| ), slower accretion 
of separate subclumps l|Larsonn^69^^'S eari^^Zim^^78jl . and dissolution of globular 
clusters ( [Aguilar, Hut, fc Ostrikcr 1988| ). More recent hierarchical models suggest that 
spheroids (bulges and halos) form in a repetitive process during the mergers of galaxies 
and protogalaxies, while disks form by slow accretion of gas between merging events (e.g., 
IS^ hitc fc Frcnk 1991). The discovery of the Sgr dwarf galaxy ( |Ibata, Gilmore, fc Irwin 1994| ) 
sparked renewed interest in halo formation through accretion of dwarf galaxies, and 
ambitious programs are now underway to map out the spatial distribution, kinemat- 
ics, and chemical abundance in the halos of the Milky Way (e.g., Morrison et al. 2000j 
[Majewski et al. 2000| ) and Andromeda (e.g., [Ferguson et al. 2002| ). In the meantime, the 
realization that hierarchical models based on cold dark matter almost inevitably predict 
many more dwarf galaxies than are actually seen around the Milky Way (,Moore et al. 1999| 
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has led to suggestions that most of the dwarf galaxies formed in the early Universe have 
dissolved into the halo (e.g., |Bullock, Kravtsov, fc Weinberg 200d| . Whether or not such 
accretion is the dominant source of stars in the stellar halo, it is likely that dwarf galaxies 
do contribute, and at large galactocentric distances their stars can remain in coherent 
orbital streams for many Gyr. Indeed, one such stream has been found in the halo of 
Andromeda Ijlbata et al. 20011 l^cConnachie et al. 2003|l : Ferguson et al. (2002) demon- 
strated dramatic evidence for large-scale debris trails and substructure in the M31 halo 
and outer disk, yet their deep imaging and star counts also found an apparently smooth 
component extending to large radii. 

The most direct way to measure ages in a stellar population is to construct a color- 
magnitude diagram (CMD) that reaches to stars below the main sequence turnoff. For 
decades, researchers have used such data to determine the ages of Galactic globular clus- 
ters and satellite galaxies of the Milky Way, but until now there has been no instrument 
capable of resolving these stars in a massive galaxy outside of our own. With the instal- 
lation of the Advanced Camera for Surveys (ACS: lR)rd et al. 1998il on the Hubble Space 
Telescope (HST), this is now possible. We have obtained deep ACS observations of a 
minor axis field in Andromeda's halo, w 51 arcmin (11 kpc) from the nucleus. In these 
proceedings, I present the preliminary results of our program IjBrown et al. 2003jl . I start 
with a brief review of the age diagnostics available in a deep CMD, then describe our 
ACS observations and their analysis, and finish with the implications for the formation 
history of Andromeda. 

2. Age Diagnostics 

Stellar population ages are best determined from resolved optical photometry reaching 
the main sequence. Although the location of the main sequence turnoff in the CMD is the 
primary age indicator, photometry reaching well below the turnoff is required to define 
the turnoff accurately. Ages among the Galactic globular clusters are usually determined 
via one of two diagnostics (figure^: the luminosity difference between the turnoff and 
the horizontal branch ( |Sandage 1982|IIben fc Renzini 1984|l . and the color difference be- 
tween the turnoff and the red giant branch ( |VandenBerg, Bolte, fc Stetson 1990| ). The 
luminosity of the horizontal branch and the color of the red giant branch are both rela- 
tively insensitive to age (at fixed composition), while the main sequence turnoff becomes 
both fainter and redder at increasing age. Our ACS observations used the F606W (broad 
V) and F814W (/) filters; for roughly every Gyr beyond 10 Gyr, the main sequence turnoff 
becomes fainter by ~ 0.1 mag in mpsi^w and redder by 0.01 mag in mpQOQw —"mFsiiw 
(see figure nj. 

The age-metallicity degeneracy is a major uncertainty when characterizing individual 
stars via photometry, or unresolved populations via low- resolution spectra (figure EJt) . 
However, the situation is better for a resolved stellar population. Because the red giant 
branch is much more sensitive to metallicity than age, and the subgiant branch and main 
sequence turnoff are sensitive to both age and metallicity, the distribution of stars on 
the red giant branch, subgiant branch, and main sequence turnoff can be simultaneously 
reproduced only for specific combinations of age and metallicity (figure EJ. 

3. Observations 

Using the ACS Wide Field Camera (WFC), we obtained deep optical images of a field 
along the southeast minor axis of the M31 halo, at 02000 = 00'*46™07*, (S2000 = 40°42'34" 
(figureEI). The field was previously imaged by Holland et al. (1996) with the Wide Field 
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Figure 1. The variation in optical photometry of a stellar population with age at fixed metal- 
licity (labeled). Isochrones at 8, 10, and 12 Gyr (labeled) are shown in a CMD constructed using 
the ACS bandpasses F606W and F814W. Note that the subgiant branch (SGB) and main se- 
quence turnoff (MSTO) are more sensitive to age than the red giant branch (RGB). The primary 
age indicators in the CMD are the luminosity difference between the MSTO and the horizontal 
branch (HE), which becomes larger as age increases, and the color difference between the MSTO 
and the RGB, which becomes smaller as age increases. However, in order to clearly define the 
MSTO, photometry must extend to the main sequence (MS) stars below the MSTO. 



Planetary Camera 2 (WFPC2); the field is not associated with the tidal streams and 
substructure found by Ferguson et al. (2002), and lies just outside the "flattened inner 
halo" in their maps. Given the nearly edge-on disk, the contribution of disk stars at this 
position should be < 3% IjWalterbos fc Kennicutt 19881 Holland et al. 1996 and references 
therein). We chose this field to optimize the crowding (trading off population statistics 
versus photometric accuracy) and to place an interesting M31 globular cluster (GC312; 
[Sarg ent et al.J j77[) near the edge o f our images. The metallicity of GC312 ([Fe/H] = 
-0.7; Huchra ^Brodie, fc Kent 1991D is near the peak in the metallicity distribution for 
the M31 halo; this should simplify our attempts to derive relative ages for this cluster 
and the halo (currently underway). Due to scheduling and orientation constraints, one 
bright foreground star [V ~ 14 mag) was unavoidable; the star and its window reflection 
affect a few percent of the total image area, which we discard. 

From 2 Dec 2002 to 11 Jan 2003, we obtained 39.1 hours of ACS images in the F606W 
filter (broad V) and 45.4 hours in the F814W filter (/), with every exposure dithered to 
allow for hot pixel removal, optimal sampling of the point spread function, smoothing of 
the spatial variations in the detector response, and filling the gap between the two halves 
of the 4096 x 4096 pix detector. We co-added the M31 images using the IRAF DRIZZLE 
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Figure 2. Panel a: An isochrone of intermediate age and metallicity {dark curve, labeled) and an 
older, more motal-poor isochrone {light curve, labeled) intersect, demonstrating that two-color 
photometry for a single star does not uniquely determine age and metallicity. Both isochronos 
also have the same color at the main sequence turnoff, and thus unresolved populations with 
these parameters would appear very similar. However, resolved photometry could distinguish 
between two such populations. Panel b: Isochrones spanning 8-12 Gyr at 1 Gyr intervals, with 
a range in metallicity {labeled) anticorrelatcd with age, such that the youngest isochrone is the 
most metal-rich, as might be expected for a simple model of chemical evolution. The subgiant 
branch and main sequence turnoff both become fainter and redder at increasing age or increasing 
metallicity, while the red giant branch is more sensitive to metallicity than age. Because the 
effects of age and metallicity counteract each other at the subgiant branch, it appears very 
narrow compared to the red giant branch. Panel c: The same metallicities as in 6, but all of 
the isochronos have the same age. Panel d: The same metallicities and ages as in b, but now 
age is correlated with metallicity, such that the older populations are more metal-rich. This is 
somewhat unpliysical, but demonstrates a situation of renewed star formation from the infall 
of primordial gas. Note the width of the subgiant branch compared to the red giant branch. 
Panels b-d demonstrate that resolved photometry of a mixed population can disentangle age 
and metallicity, because each evolutionary phase responds differently to these parameters. 
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Figure 3. A 50' x 60' Digital Sky Survey image of Andromeda, showing the positions of our 
ACS and WFPC2 images (labeled) The A CS field is 50 arcmin (11 kpc) from the center of 
M31. The nearly edge-on (i « 12.5°: I3e Vaucoulcurs 1958J disk should contribute < 3% of the 
stars in the ACS field. 

package, with masks for the cosmic rays and hot pixels, resuhing in geometricaUy-correct 
images with a plate scale of 0.03" pixel^^ and an area of approximately 210" x 220" 
(figure We then performed both aperture and PSF-fitting photometry using the 
DAOPHOT-II package IjStetson 1987|l . assuming a variable PSF constructed from the 
most isolated stars in the images. The aperture photometry on isolated stars was cor- 
rected to true apparent magnitudes using TinyTim models of the HST PSF l|Krist 1995|l 
and observations of the standard star EGGR 102 {a, V = 12.8 mag DA white dwarf) in 
the same filters, with agreement at the 1% level. The PSF-fitting photometry was then 
compared to the corrected aperture photometry, in order to derive the offset between 
the PSF-fitting photometry and true apparent magnitudes. Our photometry is in the 
STMAG system: m = — 2.5x logio/A — 21.1. For those more familiar with the Johnson 
V and Cousins / bandpasses, a 5,000 K stellar spectrum has V — mpgoew = —0.05 mag 
and / — mpsiiw = —1-28 mag. 

Of the ~300,000 stars detected, we discarded those within the GC312 tidal radius (10"; 
[Holland, Fahlman, fc Richer 1997| ), within 14.5" of a bright foreground star, within 12.6" 
of this star's window reflection, and near the image edges, leaving Ri223,000 stars in the 
final catalog. Using the SExtractor code fBertin fc Arnouts 1996;! . we estimate < 5% 



6 



T.M. Brown: The Star Formation History in the Andromeda Halo 



rN 



,^jf:v. 




f '^*"^- ■ ■'• ' •^^■'■■^/■■"^^-■'k" ^^rflHP^v!-.'.!;= ';^y;::i^">S:-L '! ' ■'/■Vs 



Figure 4. The ACS F606W image of the Andromeda halo, shown at a log stretch, and binned 
for display purposes by a factor of 20, to 0.6" pixel"^. The field subtends 211" x 221" (not much 
larger than the WFC field of view, because the largest dither was about 6"). For our preliminary 
analysis, we excluded the area around GC312 {bottom) and around a bright {V — 14 mag) fore- 
ground star {upper right) and its window reflection (to the right of the star). Many background 
galaxies can be seen through the 300,000 stars detected in Andromeda's halo, but they are a 
small (< 5%) contamination. The image has been corrected for the strong geometric distortion 
in the camera, giving the field its unusual shape. 



of the stars are contaminated by extended sources (figure [SJ. Extensive artificial star 
tests determine the photometric scatter and completeness as a function of color and 
luminosity. The CMD shows no obvious differences when comparing the population in a 
10-100" annulus around GC312 to that beyond 100"; the cluster does not appear to be 
associated with an extended underlying system. By integrating our catalog, we estimate 
that the surface brightness in our field is ~ 26.3 mag arcsec"^. 

We also obtained coordinated parallel WFPC2 observations of a second field along 
the minor axis of Andromeda (figure OJ using the same bandpasses as those in the ACS 
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Figure 5. Subsections (15" x 15") of the ACS image show the globular cluster GC312 in 
Andromeda's halo {lower right), background galaxies, and the level of crowding in the halo 
population. The images have been rebinned to a scale of 0.06" pixel"^ for display purposes. 

observations. The WFPC2 data are not nearly as deep as the ACS data, but they are 
about as deep as the Hubble Deep Field l|Willianis et al. 1996|l . and thus much deeper 
than any Andromeda observations prior to our program. Although the WFPC2 images 
resolve stars at the main sequence turnoff, the resulting CMD is not deep enough to 
characterize the turnoff and subgiant branch well. The striking differences between the 
WFPC2 and ACS CMDs are a testament to the technical advances achieved by the HST 
servicing missions (figure EJ. 

Because the ACS is a new instrument, our program includes ACS observations of 
five Galactic globular clusters spanning a wide metallicity range (Table 1), using the 
same filters utilized in our Andromeda halo observations. These cluster images allow 
the construction of empirical isochrones in the ACS bandpasses, which can be directly 
compared to the Andromeda CMD and used to calibrate the transformation of the the 
theoretical isochrones to the ACS bandpasses. Figure [3 shows the ACS F606W image 
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Figure 6. Left panel: The CMD constructed from the ACS images of the Andromeda halo, using 
aperture photometry. There are too many stars to plot them individually; instead, this is a Hess 
diagram at a logarithmic stretch. Right panel: The CMD constructed from the parallel WFPC2 
images, in a field a few arcmin further out along the minor axis of Andromeda, again using 
aperture photometry. The filters and exposure times in the WFPC2 data are approximately the 
same as those in the ACS data. Although the WFPC2 data reach stars at the main sequence 
turnoff, they are not deep enough to characterize the subgiant branch and the turnoff well. 

of M92. M92, NGC6752, and 47 Tuc are the most useful calibrators in our program, 
because their parameters are known very well; NGC5927 and NGC6528 are also useful 
because of their high metallicities, but their parameters are less secure, and they suffer 
from high, spatially variable foreground reddening (jHeitsch &: Richtler 1999|l . To increase 
the dynamic range, three images were taken in each bandpass for each cluster, with the 
exposure times varying by an order of magnitude. To minimize the number of orbits 
required for the cluster data, these observations were not dithered. Thus, we drizzled 
the data without plate scale changes, in order to remove cosmic rays and to correct 
for geometric distortion. The cluster images are significantly less crowded (on average) 
than those of M31, and the 0.05" pixels undersample the PSF, so we performed aperture 
photometry but no PSF fitting, and corrected the aperture photometry to true apparent 
magnitudes (figure |Hll. 

4. Analysis 

The ACS CMD reveals, for the first time, the main sequence population in the M31 
halo. The horizontal branch extends from a well-populated red clump to a minority blue 
population (~ 10% of the total horizontal branch population). The horizontal branch is 
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Figure 7. The ACS F606W image of M92, shown at a log stretch and binned by 20 to 0.6" 
pixel"^ for display purposes. The globular cluster images in our program were not dithered, so 
a gap appears between the two WFC chips. 



not noticeably extended - the hot horizontal branch stars that are seen in clusters like 
NGC6752 are missing. The broad red giant branch indicates a wide metallicity distri- 
bution extending to near-solar metallicities, long- known to be characteristic of the M31 
halo IjMould fc Kristian 1986|l . The luminosity function "bump" on the red giant branch 
is another metallicity indicator, becoming fainter as metallicity increases; it slopes away 
from the red horizontal branch until the luminosity difference reaches saO.5 mag - an- 
other indication of near-solar metallicities. There is also a prominent blue plume of stars 
significantly brighter than the main sequence turnoff; this minority population (~ 2% the 
size of the population at the turnoff ±1 mag) may include binaries, blue stragglers, or a 
residual young stellar population. Although the blue horizontal branch stars are charac- 
teristic of very old, metal-poor globular clusters (such as M92), the luminosity difference 
between the turnoff and horizontal branch is smaller than expected for a purely old stellar 
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Figure 8. The CMD created from our ACS images of M92. The ridge line for this cluster {dark 
grey curve) was created by taking the median mpeoew — rnpsuw color and median mpiuw 
for a series of points along the main sequence and red giant branch. A 14 Gyr isochrone at 
[Fe/H]=— 2.14 (light grey curve), transformed to the ACS bandpasses, agrees well with the ridge 
line over the CMD region used for fitting the Andromeda halo population. 

population. This is shown in figure El which shows the ridge lines and horizontal branch 
loci for the five globular clusters we observed with ACS, superimposed upon the CMD of 
Andromeda. The M31 subgiant branch is nearly horizontal, indicating a high metallicity, 
while its ridge is appreciably brighter than those of 47 Tuc, NGC5927, and NGC6528, 
implying the presence of a significantly younger population in the M31 halo. 

The comparison between the globular clusters and Andromeda in figure [HI provides a 
completely empirical indication of the age spread in Andromeda's halo. M92, NGC6752, 
and 47 Tuc were shifted in color and magnitude according to the differences in reddening 
and distance between the clusters and Andromeda. Those shifts naturally aligned the 
horizontal branch loci of these clusters to the horizontal branch of Andromeda, thus 
demonstrating the accuracy of the parameters in Table 1, but we could have simply 
aligned the horizontal branches without any knowledge of the relative distances and 
reddenings. Indeed, forcing alignment at the horizontal branch was required for NGC5927 
and NGC6528, because the parameters of those clusters are very uncertain. Once the 
clusters are all aligned at the horizontal branch, the red giant branches of the clusters 
span the broad red giant branch of Andromeda, empirically demonstrating the wide 
metallicity distribution in Andromeda's halo. For the most metal-poor cluster (M92; 
lightest grey curve in figure El, the subgiant branch luminosity agrees well with that in 
Andromeda, and the turnoff of the cluster agrees well with the blue edge of the turnoff in 
Andromeda; this indicates that Andromeda contains a significant population of metal- 
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Figure 9. Left panel: The CMD of Andromeda constructed from the PSF-fitting photome- 
try of the ACS images. Completeness hmits determined from artificial star tests are marked. 
Right panel: The CMD of Andromeda, with the ridge lines {grey curves) and horizontal branch 
loci {grey diamonds) of five Galactic globular clusters superimposed. The shading of the ridge 
lines and horizontal branch points range from light grey to dark grey as the cluster metallicity 
increases (see Table 1). The data for M92, NGC6752, and 47 Tuc have only been shifted by 
the differences in distance and reddening between the clusters and M31, yet their horizontal 
branch loci agree well with the horizontal branch of M31. Those parameters are very uncertain 
for NGC6528 and NGC5927, so the data for the two metal-rich clusters were shifted to align 
their horizontal branch loci with the M31 horizontal branch: NGC6528 was shifted 0.16 mag 
brighter, while NGC5927 was shifted 0.11 mag brighter and 0.05 mag redder. Moving from the 
most metal-poor cluster (M92) to the most metal-rich (NGC6528), the red giant branches of 
the clusters span the width of the red giant branch in M31, yet the subgiant branches of the 
clusters become increasingly fainter than that of M31, indicating that the metal-rich stars in 
the M31 halo are much younger than those in the clusters. 



poor stars with ages similar to the oldest Galactic globular clusters. However, as one 
compares clusters at increasing metallicity to Andromeda, discrepancies at the subgiant 
branch and turnoff become more apparent. Although the red giant branches of 47 Tuc, 
NGC5927, and NGC6528 straddle a large fraction of the red giant branch population 
in Andromeda, the cluster subgiant branches are well below the bulk of the subgiant 
population in Andromeda. This is a strong indication that the metal-rich populations in 
Andromeda are significantly younger than the cluster, given that the subgiant branch 
and turnoff become fainter by approximately 0.1 mag for every 1 Gyr increase in age. 

Our full analysis of the star formation history is in progress. Here, I will supplement 
the comparisons shown in figure El with examples from the modeling to date (figure 
Our modehng is based upon the isochrones of VandenBerg, Bergbusch, & Dowler (2003, 
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Table 1. M31 and globular cluster parameters 



Exposure 



Name 


(m - M)v 


Eb-v 






) 


Name 


(mag) 


(mag) 


[Fe/H] 


F606W 


F814W 


M31 


24.68'' 


0.08*= 


-0.6^= 


140870 


163570 


M92 




0.023*^ 


-2.14'' 


95.5 


106.5 


NGC6752 


13.17'= 


0.055*= 


-1.54^ 


44.5 


49.5 


47 Tuc 


13.27*5 


0.032*= 


-0.83' 


76.5 


78 


NGC5927 


15.81'' 


0.45'' 


-0.37*' 


532 


355.7 


NGC6528 


16.15' 


0.55' 


-0.2' 


504 


371 



''Freedman & Madore (1990). ""Schlegel et al. (1998). '=Mould & Kristian (1986). ''VandenBerg & 
Clem (2003). ''Renzini et al. (1996). 'VandenBerg (2000). ^^Zoccali et al. (2001). ''Harris (1996). 
'Momany et al. (2003). 



in preparation) , which show good agreement with cluster CMDs spanning a wide range of 
metallicity and age (e.g., figure IS)). To transform these isochrones from the ground-based 
bandpasses to the HST-based bandpasses, we used the spectra of Lejeune, Cuisinier, & 
Buser (1997) to calculate V — mpQi^Qw and I — mp&iAW , as a function of effective tempera- 
ture and gravity along the isochrones, then applied those differences to the ground-based 
magnitudes of the isochrones, with a smaU (< 0.05 mag) empirical color correction to 
force agreement with our globular cluster CMDs. After this correction, the isochrones 
match the ridge lines within < 0.02 mag over the region of the CMD we are fitting. 
Our observed CMDs of these clusters are reproduced by a 12.5 Gyr isochrone for 47 Tuc 
and 14 Gyr isochrones for NGC6752 and M92. The isochrones do not include core He 
diffusion, which would reduce their ages by 10-12%, thus avoiding discrepancies with the 
age of the Universe ( [VandenBerg et al. 2002| ) . 

Using isochrones with a range of ages and metallicities, we fit the region of the M31 
CMD shown in figure [TUk using the StarFish code of Harris & Zaritsky (2001). Restrict- 
ing the fit to this region of the CMD focuses on the most sensitive age and metallicity 
indicators while avoiding regions of the CMD that are seriously incomplete, sparsely 
populated, or poorly constrained by the models (e.g., the horizontal branch morphol- 
ogy). Note that we do not include the red giant branch luminosity function bump in our 
fitting; although theory predicts that this bump becomes fainter with increasing metal- 
licity, the theoretical zeropoint is fairly uncertain. The StarFish code fits the observed 
CMD through a linear combination of input isochrones, using minimization, where 
the isochrones are scattered according to the results of the artificial star tests. We first 
followed the standard method used when researchers determine the metallicity distribu- 
tion from the red giant branch (figure ITUb ). We fit the red giant branch using a set of 
old (13.5 Gyr) isochrones spanning a wide metallicity range (—2.31 < [Fe/H] < 0). The 
resulting metallicity distribution was similar to that found by Holland et al. (1996) in 
our same field, and Durrell et al. (2001) in a field 20 kpc from the nucleus. However, it 
is very obvious that these purely old isochrones do not match the subgiant branch or 
the turnoff. Next, as shown in figure 11111 :. we tried fitting the entire region of figure [Tnb , 
with a wider age range (11.5-13.5 Gyr), but found no acceptable combination. The insets 
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Figure 10. Panel a: The region of the Andromeda halo CMD that we used for fitting the star 
formation history. Restricting the fit to this region avoids parts of the CMD that are seriously 
incomplete, sparsely populated, or poorly constrained by the theoretical models. A histogram 
(inset) of the number of stars along a cut through the subgiant branch {grey line) highlights the 
differences between the data and the models in the subsequent panels. Panel b: The fit to the red 
giant branch using only 13.5 Gyr isochrones spanning a wide range in metallicity. Although such 
models can reproduce the width of the Andromeda red giant branch, the subgiant branch and 
turnoff are much fainter than those in the data. The histogram (inset) compares the number of 
stars along the cut (black line) through the subgiant branch to that in the data (grey histogram); 
they are poorly matched. Panel c: The fit to the entire region shown in panel a, using isochrones 
with a wide metallicity range but ages of 11.5-13.5 Gyr. The histogram (inset) again highlights 
that this model does not reproduce the distribution seen in the data. Panel d: The fit to the data 
shown in panel a, using a wide metallicity range and ages of 6-13.5 Gyr. This model reproduces 
the data well, as highlighted by the histograms (inset). 

in figure 1101 contain histograms showing the number of stars (data: grey] model: black) 
along a cut through the subgiant branch (thick line). In figure [TUt . the residual subgiant 
stars that are not reproduced by this old model (inset) suggest that at least 20% of the 
stars in the halo must be younger than 11.5 Gyr. We conclude that a purely old stellar 
population cannot explain the CMD of the M31 halo. 
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Age (Gyr) 

Figure 11. Upper panel: The star formation history of the metal-rich stars in the best model 
(figure llOH 'l obtained to date. The bulk of the metal-rich population is of intermediate age. 
Lower panel: The star formation history of the metal-poor stars in the best-fit model. Most of 
the metal-poor stars are old, but this model includes a small population of young metal-poor 
stars. Note that this star formation history is not unique; we can obtain a fit that is nearly as 
good as this fit, using two completely distinct population components (a purely old metal-poor 
population and a purely young metal-rich population). The detailed star formation history will 
be explored more fully in a future paper. 

Next, we expanded the age range to 6-13.5 Gyr and repeated the fit (see figure ITUli') . 
The width of the red giant branch is now matched without a mismatch at the subgiant 
branch (inset). The best-fit model can be broadly characterized by a combination of two 
dominant populations (figure [TT|l : 56% of the stars (red) are metal-rich ([Fe/H]> —0.5) 
and of intermediate age (6-11 Gyr), while 30% of the stars (blue) are metal-poor ([Fe/H]< 
—0.5) and old (11-13.5 Gyr). About half of this metal-rich population (i.e., 28% of the 
total population) is 6-8 Gyr old. Note that these models only illustrate, in a broad sense, 
the dominant populations present in the M31 halo. Other combinations of young metal- 
rich and old metal-poor stars are possible. For example, we produced a similar fit by 
combining two very distinct isochrone groups: one at 6-8 Gyr with [Fe/H]> —1, and one 
at 11.5-13.5 Gyr with [Fe/H]< —1. We are currently investigating detailed constraints 
on the star formation history. 

5. Disk Contamination 

Note that our halo field should be relatively free of stars moving in the stellar disk. 
Several estimates for the disk contribution at this distance range from 1-3%, even with 
a significant thick disk component (Holland et al. 1996 and references therein). On aver- 
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Figure 12. Top panel: A simple schematic showing the intersection of our sight-line with 
the Andromeda disk. The plane of the disk is inclined 12.5° with respect to our line of sight 
{dc Vaucoulcurs 1958), so that we are viewing the halo population at 11 kpc from the nucleus 
but the disk population at 51 kpc from the nucleus. The disk has an exponential scale length of 
5.3 kpc (Waltcrbos & Kcrmicutt 19881. Note that the southeastern half of the disk is the side 
furthest from the observer (lyc & Richtcr 19851. Middle panel: The isophotes for Andr omeda 
show evidence of a warp near 22 kpc on the major axis (|Walterbos fc Kennicutt 1987^ . Here, 
I insert a hypothetical 10° warp at 22 kpc on the minor axis, such that more of the disk pop- 
ulation is moved into our sight-line. Our sight-line has effectively moved inward by about 2 
scalelengths, increasing the disk population by a factor of 8. However, even with this large warp 
into our sight-line, the increase in the disk population would not explain the large population of 
intermediate-age stars we find in our field. Bottom panel: A more realistic depiction of the warp 
in the Andromeda disk, by interpolating the measurements of Braun (1991) along the minor 
axis. Early models of the warp (e.g., Henderson 19791 imply our sight-line would not intersect 
the disk at a ll, because of the increasing inclination at large distances from the nucleus. Later 
models (e.g.. f5raun 19911 indicate that the inclination varies from 25° close to the nucleus to 
10-15° further out, as depicted here. The result is that our effective position within the disk is 
significantly beyond the 51 kpc shown in the top panel. 



age, the plane of the disk is inclined from our sight-line by 12.5° l|de Vancouleurs 19581 
figure El top panel), so we are observing halo stars 11 kpc from the nucleus but disk 
stars 51 kpc from the nucleus. A warp in the outer disk, if it angled the disk into our 
sight-line, would move the sight-line intersection with the disk inward, but even a large 
warp could not explain the number of young stars seen in our field. E.g., Walterbos & 
Kennicutt (1987) found a possible warping of the stellar disk at 100' (22 kpc) on the 
southwest major axis. If we assume that the stellar disk along the southeast minor axis 
is warped such that the disk beyond 22 kpc is inclined by an additional 10° into our 
sight-line, this would move our effective disk position inward by about two scalelengths 
(figure [T21 middle panel), increasing the disk population by a factor of 8 (i.e., the disk 
would contribute wlO-20% of the total population in our field). However, we find that 
56% of the stars in our field are metal-rich and of intermediate age, so this hypothetical 
increase in the disk contamination would not account for our surprisingly young popula- 
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tion. In reality, studies of the neutral gas kinematics fe.g.. lBraun 199lllHenderson 1979)l 
find that the outer disk is actually viewed significantly closer to edge-on than the inner 
disk f figure IT^ bottom panel), such that the warp in the disk would move our sight-line 
further out, beyond the intersection given by a flat disk. Finally, the fact that the metal- 
licity distribution in our field is very similar to that twice as far out (Holland et al. 1996; 
[Durrell, Harris, &: Pritchet 2001| ) provides another strong indication that we are viewing 
a halo-dominated population. 



6. Summary and Discussion 

The CMD of the M31 halo is evidently inconsistent with a population composed solely 
of old (globular cluster age) stars; instead, it is dominated by a population of metal-rich 
intermediate-age stars. Although the high metallicity in the M31 halo is well-documented, 
the large age spread required to simultaneously reproduce the red giant branch, subgiant 
branch, and main sequence came to us as a surprise. Earlier studies of the red giant branch 
were insensitive to this age spread. For example, Durrell et al. (2001) were able to explain 
the metallicity distribution 20 kpc from the nucleus, with a simple chemical evolution 
model forming most of the stars at very early times (see also lCote et al. 2000|l . Although 
our field is relatively small in sky coverage, it appears representative; the metallicity in 
our field (Holland et al. 1996) agrees well with that much further out (Durrell et al. 2001), 
and there are no indications of substructure or tidal streams in the region we surveyed 
HFerguson et al. 2002) . 

It seems unlikely that star formation in the halo proceeded for ~6 Gyr from gas in 
situ; instead, the broad age dispersion in the halo is likely due to contamination from 
the disruption of satellites or of disk material into the halo during mergers. Indeed, our 
current analysis of the data cannot rule out the possibility that M31 and a nearly-equal- 
mass companion galaxy experienced a violent merger when the Universe was half its 
present age. If the 6-8 Gyr population in the halo represents the remnants of a disrupted 
satellite, the relatively high metallicity suggests that it must have been fairly massive. On 
the other hand, the stars could represent disruption of the M31 disk, either by a major 
collision when M31 was ~ 6 Gyr old, or by repeated encounters with smaller satellites. 
The resulting halo would be a mix of the old metal-poor stars formed earliest in M31's 
halo, disk stars that formed prior to the merger(s) that were subsequently dispersed 
into the halo, stars formed during the merger (s), and the remnant populations of the 
disrupted satellite(s). 

This work was done in collaboration with H.C. Ferguson, E. Smith (STScI), R.A. Kim- 
ble, A.V. Sweigart (NASA/GSFC), R.M. Rich, D. Reitzel (UCLA), A. Renzini (ESQ), 
and D.A. VandenBerg (U. of Victoria). I am grateful to J. Harris (STScI) and P. Stetson 
(DAO) for respectively providing the StarFish and DAOPHOT-II codes and assistance 
in their use. Thanks to the members of the scheduling and operations teams at STScI 
(especially P. Royle, D. Taylor, and D. Soderblom) for their efforts in executing a large 
program during a busy HST cycle. Support for program 9453 was provided by NASA 
through a grant from the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 
5-26555. The Digitized Sky Survey was produced at the Space Telescope Science Institute 
under U.S. Government grant NAG W-2166. The images of these surveys are based on 
photographic data obtained using the Oschin Schmidt Telescope on Palomar Mountain 
and the UK Schmidt Telescope. The plates were processed into the present compressed 
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